Newtonian polytropes for anisotropic matter: General framework and applications 
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We set up the general formalism to model polytropic Newtonian stars with anisotropic pressure. 
We obtain the corresponding Lane-Emden equation. A heuristic model based on an ansatz to 
obtain anisotropic matter solutions from known solutions for isotropic matter is adopted to illustrate 
the effects of the pressure anisotropy on the structure of the star. In particular, we calculate 
the Chandrasekhar mass for a white dwarf. It is clearly displayed how the Chandrasekhar mass 
limit changes depending on the anisotropy. Prospective astrophysical applications of the proposed 
approach are discussed. 
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I. INTRODUCTION 

In the context of Newtonian gravity, polytropic equa- 
tions of state are particularly useful to describe a great 
variety of situations (see Refs. [TUB] and references 
therein), their great success stemming mainly from the 
simplicity of the equation of state and the ensuing main 
equation (Lane-Emden). Polytropes in the context of 
general relativity have been considered in Refs. [71113] 
(and references therein). However, in this work, we re- 
strict the analysis to Newtonian polytropes. 

The theory of polytropes is based on the polytropic 
equation of state 



P = Kp'^ = Kp^+^''- 



(1) 



where P and p denote the isotropic pressure and the mass 
(baryonic) density, respectively. Constants K, 7, and 
n are usually called the polytropic constant, polytropic 
exponent, and polytropic index, respectively. 

The polytropic equation of state may be used to model 
two very different types of situations, namely: 

(i) When the polytropic constant K is fixed and can 
be calculated from natural constants. This is the 
case of a completely degenerate gas in the nonrel- 
ativistic (7 = 5/3; n = 3/2) and relativistic limit 
(7 = 4/3; n = 3). Polytropes of this kind are par- 
ticularly useful to model compact objects such as 
withe dwarfs (WDs), and they lead in a rather sim- 
ple way to the Chandrasekhar mass limit. 

(ii) When K in a, free parameter, as, for example, in 
the case of isothermal ideal gas, or in a completely 
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convective star. Models related to isothermal 
ideal gas are relevant in the so-called Schonberg- 
Chandrasekhar limit (see Ref. [4] for details). 

Our motivation to extend polytropic stellar models to 
cases in which the pressure anisotropy is allowed is based 
on the fact that the local anisotropy of pressure may be 
caused by a large variety of physical phenomena of the 
kind we expect in compact objects (see Ref. [16] for an 
extensive discussion on this point). 

Indeed, the study of anisotropic (principal stresses un- 
equal) spherically symmetric fluids has a long and ven- 
erable story. It started with Jeans [T7], who studied the 
anisotropy produced by anisotropic velocity distributions 
in galaxies. The first mention of local anisotropy of pres- 
sure in spherically symmetric selfgravitating fluids may 
be found in the seminal paper by Lemaitre [18] . In page 
63 of that paper, Lemaitre realizes that the stringent 
limit in the compactness of a homogeneous relativistic 
sphere is related to the isotropy of pressure, and there- 
fore he proposes to relax that condition. He considers 
the "limiting" case in which the radial pressure vanishes, 
but the tangential does not. However, the interest in this 
subject started to grow exponentially after the pioneering 
work of Bowers and Liang 19'. For recent references on 
this subject, see Refs. [20-24 (and references therein). 
An alternative approach to anisotropy comes from kinetic 
theory using the spherically symmetric Einstein- Vlasov 
equations, which admits a very rich class of static solu- 
tions, none of them isotropic (Ref. [25] (and references 
therein). The advantages or disadvantages of either ap- 
proach are related to the specific problem under consider- 
ation. As we shall see below, our method links our models 
continually with the isotropic case (see Sec. IV), thereby 
allowing us to bring out the influence of anisotropy on 
the structure of the object. Evidently, both approaches 
should give the same physical results. 

Among all possible sources of anisotropy, there is one 
particularly related to our endeavor in this manuscript. 



namely the intense magnetic field observed in compact 
objects such as white dwarfs, neutron stars, or magne- 
tized strange quark stars (see, for example, Refs. [^BHSO] 
and references therein). 

Indeed, it is a well-established fact that a magnetic 
field acting on a Fermi gas produces pressure anisotropy 
(see Refs. [5TH55] and references therein). In some way, 
the magnetic field can be addressed as a fluid anisotropy. 

Particularly appealing is the fact that magnetic fields 
may severely affect the Chandrasekhar mass limit of a 
white dwarf [36 38] . 

For all the reasons above, we intend in this paper to 
develop the general formalism to describe polytropes in 
the presence of pressure anisotropy. 

For the sake of completeness, we shall first review very 
briefly the theory of polytropes for a perfect (isotropic) 
fluid. Next, we shall display the general formalism for 
anisotropic fluids. In order to bring out the effects of 
anisotropy on the structure of the star, we shall further 
assume an ansatz allowing us to calculate the influence 
of anisotropy on the Chandraskhar mass limit. Finally, 
we shall conclude with some possible applications and 
unanswered issues. 



Pc 



(7) 



where the subscript c indicates that the quantity is evalu- 
ated at the center, and the following boundary conditions 
apply: 



dz 



(z = 0) = 0; uj{z = 0) = 1. 



The boundary surface of the sphere is defined by z = z„, 
such that uj{zn) — 0. As is well known, bounded configu- 
rations exist only for n < 5, and analytical solutions may 
be found for n = 0, 1 and 5. 

In the case 7 = 1, which corresponds to an isothermal 
ideal gas, the ensuing Lane-Emden equation reads 



a uj 
~d^ 



2duj 

z dz 



where now 



A' 



(8) 



(9) 



II. THE POLYTROPE FOR FLUIDS WITH 
ISOTROPIC PRESSURE 



A' 



AnGpc 
K 



(10) 



Polytropes are assumed to be in hydrostatic equilib- 
rium (for deviations from this condition see Refs. [3| and 
[5]); therefore, the two starting equations are the equa- 
tion of hydrostatic equilibrium. 



dP 

dr 



dr 



(2) 



and the Poisson equation (in spherical coordinates). 



w = — 
K 

and the following boundary conditions apply: 
duj 



(11) 



dz 



(z = 0) = 0; a;(z = 0) = 0. 



We shall next generalize the scheme above to the case 
when the pressure is no longer isotropic. 



1 d , nd<^, 

(r — 

r^ dr dr ' 



AnGp, 



(3) 



with $ and G denoting the Newtonian gravitational po- 
tential and the gravitational constant respectively. 

Combining the two equations above with Eq. (II]), one 
obtains after some simple calculations the well-known 
Lane-Emden equation (for 7 7^ 1): 



dzJ 



2duj 

z dz 



0, 



with 



z 
I' 



A' 



inGp. 



(n-l)/n 



K{n + 1) 



(4) 



(5) 



(6) 



III. THE POLYTROPE FOR ANISOTROPIC 
FLUIDS 

If we allow the principal stresses to be unequal, then 
the hydrostatic equilibrium equation reads 



dPr 
dr 



d<^ 
dr 



P 



fPl -Pr 



(12) 



where Pr and P± denote the radial and tangential pres- 
sures, respectively. This is the Newtonian limit of 
the generalized Tolman-Opphenheimer-Volkoff equation 
for anisotropic matter. We recall that the Tolman- 
Opphenheimer-Volkoff equation is a consequence of Ein- 
stein equations and Bianchi identities, or equivalently, it 
comes directly from the vanishing of the covariant diver- 
gence of the energy momentum tensor. If spherical sym- 
metry is assumed, then necessarily the nonradial stresses 
are equal: Pg = Pa, = P±, the only freedom being in 



this case that Pr ^ P± ■ Indeed, spherical symmetry sup- 
poses enough freedom to rotate Cartesian axes in order 
to guarantee P^ = Py = Pj_ and Pj — Pr 7^ P±- Of 
course, if one does not assume spherical symmetry, then 
in principle all three main stresses may be different. 

For the Poisson equation, of course, we get the same 
expression [Ref. p], as both are under spherical symme- 
try. 

We shall next assume a polytropic equation [Eq. ([l])] 
for the radial pressure Pr- Then, using Eqs. ( [l2| and 
3», we may write 



® 



dr 



-jKp 



,7-2 



dp 
dr 



2A 

r p ' 



(13) 



where A = P±^ — Pr- 

For the case 7 ^ 1, we can formally integrate the equa- 
tion above between any interior r and the surface of ra- 
dius r = r^ = constant, which gives us 



$ = F(r)-ii:(n + l)pi/", 



which can be written as 



-K{n+l) 



where 



r A 

F{r) = 2 —dr. 
Introducing the variables 



ar, 



' '"^ ^-($c-Fc)r-^ 



(n + 1)"K" 



and 



p \ '/" $ - F 



<^n-Fr: 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 



For the isothermal case which corresponds to 7 = 1, 
Eq. ( 13 1 becomes 



d$ ,^ .dp 2A 

— = -Xp-i-f + — , 
dr dr rp 



which after integration yields 



p^pJF-^)/^, 



(23) 



(24) 



where the potential was set to zero at r = and F is 
defined by Eq. ( 16 1 (with the inferior limit set to r = 0, 
of course). 

Then the corresponding Lane-Emden equation be- 
comes 



d^X 2 dX 



dz 



dz^ 



with 



ar, 



47rGpc 
K 



X 



F 
K 



K' 



(25) 



(26) 



(27) 



(28) 



It is obvious that in order to proceed further with the 
modeling of the compact object [i.e., in order to inte- 
grate Eqs. (20) or (25)], we need to prescribe the spe- 



cific anisotropy of the problem (A). Such information, 
of course, depends on the specific physical problem un- 
der consideration. Here we shall not follow that direc- 
tion; instead, we shall assume an ansatz already used 
in the modeling of relativistic anisotropic stars [39l |40] , 
whose main virtue (besides its simplicity) is the fact that 
the obtained models are continuously connected with the 
isotropic case. 



IV. MODELING ANISOTROPIC POLYTROPES 



where as before, the subscript c indicates that the quan- 
tity is evaluated at r = 0, the extended Lane-Emden 
equation can be written as 

where 



In order to obtain specific models, we shall here adopt 
the nonrelativistic version of the heuristic procedure 
used in Ref. [35], which allows to obtain solutions for 
anisotropic matter from known solutions for isotropic 
matter; that is, 



Cf{r)pr 



N 



(29) 



and 



X = w + Y 



Y 



F 



•i'c-F, 



(21) 



(22) 



where C is a parameter which measures the anisotropy; 
the function / and the number N are to be specific for 
each model. Following that procedure the ansatz 



fr 



N-l 



dr ' 



(30) 



0.9 






0.8 






0.7 






0.6 


n\ 




3 0.5 


\\ \ 




0.4 






0.3 


\\ 




0.2 


\ \ 


\ \ 


0.1 




\ \ \ 

''■■'■■ \ 













X 


5 










/ x^ 


4 








/' / 


/' 


3 3 




// 


/ 


/ 

/ 
/ 


/ 


2 




/ 


/ 






1 


/ / / 











0123456789 10 



10 12 14 16 



FIG. 1: li) as a function of z for n = 3 (WD) and h = 0.5 
(dashed line, upper); h = 1.0 (solid line); h = 1.5 (short- 
dashed line, lower). 



FIG. 3: w as a function of z for 7=1 (isothermal gas) and 
ft = 0.5 (dashed line, lower); ft = 1.0 (solid line); ft = 1.5 
(dot-dashed line, upper). 
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FIG. 2: Mass ratio M/Mch as a function of ft. 



leads to 



F = 2C$. 



Thus, Eq. (20 1 can be easily reduced to 



(f 



w 



2dw 

z dz 



hw'' = 0, 



(31) 



(32) 



where h — 1 — 2C. For simplicity, we assume h to be 
constant throughout the sphere, which of course does not 
imply the constancy of either pressure. Observe that 
Eq. ( 32 1 is the same as the Fowler equation [5T] when 
/3o = —h and a; = in the notation of [42] [see Eq. (2.2) 
in this last reference]. However, it should be noticed 
that Eq. (20) is more general than the Fowler- Emden 
equation. 



Now, we proceed to integrate numerically with the 
boundary conditions 

w{0)^l, J(0)-0, 

with n — 3, which represents a relativistic WD. 

Figure 1 displays the dimensionless variable w as a 
function of the dimensionless variable z for different val- 
ues of h. Figure 2 shows the Chandrasekhar mass ratio 
(with respect to the isotropic mass Mch) 



z| {dw/dz) 



M 



Mch [zj (dw/dz) 



(33) 



Z3ih=l 



as a function of the anisotropy parameter h. 

Following the same ansatz for the heuristic model in 
the case 7 = 1, the Lane-Emden equation becomes 



aw 
Iz^ 



2dw 

z dz 



^he 



where 



$ 



K 



and from Eq. ( 24 1 we get 



PcS 



-w, 



-h<S>/K 



(34) 



(35) 



(36) 



Equation (|34|) has to be integrated with the following 

dw , 



central conditions: 

w;(0) 



0, ^(0) 
dz 



0. 



Figure 3 displays the solution for different values oi h. As 
for polytropes with n > 5, the isothermal sphere consists 
of an ideal gas which has infinite radius. 

It is remarkable, although not general, to observe the 
invariance of scale of Eqs. (32) and (34 1. In fact, if we 



redefine the dimcnsionlcss spatial variable 

( — Vhz, 



(37) 



we recover the "isotropic" Lane-Emden equation, with 
all the special cases or analytic solutions for n = 0, 1, 5. 
Doing so, it is easy to see where the surface is for any h 
from calculation for h = 1. This became obvious when 
we considered the asymptotic expansion near z = (for 
7 7^ 1), rendering 



w 



1 



«'+iio«' + 



(38) 



Thus, calculations in terms of C lead to an invariant Men- 
But physically, the Chadrasekhar mass can clearly be 
stretched (or shrunk) with the anisotropy h. The same 
rescaling [Eq. (37)] works for Eq. (34) with a different 
expansion near C = 0. 



CONCLUSIONS 



of local anisotropy in such an important problem as the 
Chandrasekhar mass limit. We do not know if the in- 
ferred super-Chandrasekhar white dwarfs from collected 
data [^5^H5] are the result of anisotropy as considered 
here. But this interesting matter and the physical condi- 
tions for real stars deserve more attention elsewhere. 

We have also shown (under the same ansatz) how the 
7 = 1 case is affected by anisotropy; this might be of in- 
terest in the discussion of the Schonberg-Chandrasekhar 
limit. 

Of course, for the modeling of specific astrophysical 
objects, full information about the anisotropy (A) has to 
be provided. 

When the approach to anisotropy comes from kinetic 
theory by means of the Vlasov-Poisson system, other 
equations such as Emden-Fowler can be derived from the 
generalized polytrope equation (see Ref. 0^ and refer- 
ences therein). All this is within the context of Newto- 
nian gravity and spherical symmetry. 

Finally, we want to stress that all we have done here 
requires spherical symmetry, at least as an approxima- 
tion. It is possible that this symmetry can be broken 
by a strong magnetic field, rendering the distribution 
anisotropic and nonspherical j31j . In such a case, of 
course, the method presented here does not apply. 



We have established the general framework foe model- 
ing polytropes in the presence of anisotropic pressure. 
As mentioned in the Introduction, we undertook this 
task motivated by the conspicuous presence of such an 
anisotropy in compact objects and its influence on their 
structure. 

We also obtained some specific models based on a 
heuristic ansatz used many years ago to handle the 
anisotropy of pressure. The main purpose of that mod- 
eling was to bring out, in an explicit way, the influence 
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